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Considering that gravitational force might deviate from Newton's inverse-square law 
and become much stronger in small scale, we present a method to detect the possible 
existence of extra dimensions in the ADD model. By making use of an effective variational 
wave function, we obtain the nonrelativistic ground energy of a helium atom and its 
isoelectronic sequence. Based on these results, we calculate gravity correction of the 
ADD model. Our calculation may provide a rough estimation about the magnitude of 
the corresponding frequencies which could be measured in later experiments. 
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There are at least two seemly functional energy scales in nature, the electroweak 
scale rriEW ~ 10 2 Gev and the Planck scale M p i ~ 10 19 Gev, where gravity becomes 
as strong as the gauge interactions. Over the last two decades, explaining the small- 
ness and radiative stability of the hierarchy ijiew/M p i ~ 10~ 17 has been one of 
the greatest driving forces behind the construction of theories beyond the Standard 
Model (SM). N. Ar kani-H amed, S. Dimopoulos and G. Dvali proposed an exciting 
model (ADD model) In the ADD model, a very simple idea is to suppose there 

are n extra compact spatial dimensions of radius R and to assume that M p i^ 4+n ^ 
is around the scale of tuew- The weakness of four dimensional gravitational force 
is explained as that the origin strong high dimensional one leaks into extra dimen- 
sions. If two test massed of mass mi, m,2 placed within a distance r <C R, they will 
fell a gravitational potential dictated by Gauss's law in (4+n) dimensions 

v(r) JfH* ( r «R) (1) 

On the other hand, if the masses are placed at distances r>fl, their gravitational 
flux lines can not continue to penetrate in the extra dimensions, and the usual 1/r 
potential is obtained 

V(r)~ ™ , (r»R) (2) 
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in which the radius of extra dimensions R is determined by M p i^ +n \ as 



(3) 



If M p /( 4+Jl ) ~ ITev, one can obtain R ~ il0 _17+ ^cm. For n — 1, i? ~ 10 12 cm 
implying deviations from Newtonian gravity over solar system distances, so this case 
is empirically excluded. Particularly, the n = 2 case implies sub-millimeter extra 
dimensions, while the experimental conditions for testing Newton's inverse-square 
law (ISL) by torsion pendulum method was just about to be available when ADD's 
proposal appeared. In fact, for all n > 2, however, the modification of gravity only 
becomes noticeable at distances smaller that those currently probed by experiments 
will be performed in the very near future. Therefore, many people have devoted to 
the s earch of deviation from ISL as well as extra dimensions during the past few 
years SIS. Luo and Liu I^IZI first proposed the idea of detecting the deviation of 
Newtons inverse-square law as well as exploring extra dimensions by spectroscopy 
experiments and took hydrogen-like and muonic atomic system as illustrations. Here 
we wonder if one can detect the extra dimension by spectroscopic experiments due 
to a helium or helium-like ions. 

From the early days of the quantum mechanics the ground state ionization en- 
ergy of a helium atom was a benchmark for approximate methods of solving nonrel- 
ativistic Schrodinger equation for a few-body system. One of the earliest variational 
calculations has been performed by Hylleraas^ in 1929. In 1957, Kinoshita © ob- 
tained higher order corrections including the Lamb shift calculations confirmed a 
very good agreement with the best experimental value. We would like to mention 
here the two most recent calculations. The first is aimed to elaborate an efficient 
variational method for the many electron atoms. The second is to find an effective 
and economical way for studying a helium and helium-like two electron atoms. Sev- 
eral techniques fo r ex ample the finite-element method and the hyperspherical 
harmonic method have been used for such system, but variational method is the 
most powerful tool for studying this problem and much work has been carried out 



In this paper, we adopt a simple effective trial wave function and calculate the 
variational ground state energy of a helium atom and helium-like ions. Based on 
this trial wave function and making use of the perturbation method, we obtain 
the gravity corrections of the ADD model and the contribution to spectroscopy. 
In these calculations, we use the Mathematic software, which is a tool of symbolic 
calculation. All the results are exact except for that of the wave function cutoff, 
which would bring an error but the error can be controlled by setting the precision. 

It is known that the non-relativistic Schrodinger equation for two-electron 
helium-like systems with nuclear charge Z can be written as ( in atomic units) 



112113114115116] 



Hip(r 1 ,r 2 ) = £ty(ri,r 2 ), 



(4) 



(5) 
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where P12 are the momenta of the electrons, P = — pi — P2 is the momentum 
of the nucleus, ri and r 2 denote the positions of the electrons with respect to the 
nucleus, and the finite nucleus-to-electron mass ratio is given by M = m nuc i/me. To 
solve the above Schrodinger equation for the ground states, we employ variational 
approach by finding the stationary solutions of the following energy functional 

in which the volume element dr = iir 2 rir 2 dridr 2 dri 2 . In order to perform varia- 
tional calculation, we need to make a judicious choice for the wave function. To this 
end we adopt a simply effective variational wave function, which contains a flexible 
scaling parameter k 

^(ks,kt,ku) = e- ks Y,Ci mn {ks) l (kt) 2m (ku) n , (7) 
s = ri+?- 2 , t = — ri+T2, u = ri2, (8) 

where s, t and u are called the Hylleraas coordinates^ and the number of the terms 
of the polynomial in ip is denoted as N. Under the Hylleraas coordinates, the kinetic 
energy K, the potential energy P and the denominator W — J ijj*ipdT in Eq. ([5]) 
can be expressed as follows, respectively 

K= ds du dtTr 2 {u(s 2 -t 2 )[(d s ^) 2 

JO JO J-u 

+{d u ipf + {d t ^) 2 ] - 2 S (t 2 - u 2 )d s ^d u ^ 

+2t(s 2 - u 2 )d t i)d u ^ + 2u[{s 2 - u 2 ){d s ^) 2 
~(t 2 -u 2 )(d t ^) 2 ]/M}, (9) 

P = - ds du dtTr 2 (4Zsu- s 2 +t 2 )if> 2 , (10) 

JO Jo J-u 

W = ds du dtn 2 u{s 2 - t 2 )tj} 2 . (11) 

JO Jo J-u 

Substituting Eqs. 0, (O into Eq. ^ and using the variational method, we cal- 
culate the ground state energy of a helium atom and helium-like ions, and obtain 
the approximate wave function tp. For definition and simplicity, in the following 
calculation, we take N = 10. 

The above work is performed by a program, which is written by Mathematic 
Language. In what follows, based on the above calculating results, we investigate 
the spectroscopy correction from the ADD gravitational potential energy. Even if 
gravity would become stronger in small scale, it is still very weak compared to 
electromagnetic force. So it is convenient to treat the gravitational potential as a 
perturbation to calculate the energy correction of a helium atom and helium-like 
ions. The correction of the ground state energy is written as 

AE = 8n 2 I ^*V{r ll r 2 ,r 1 2)^r 1 r 2 dr 1 dr 2 dr 12 . (12) 
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The gravitational potential is given by 



El 



E2 G (4+ „ ) m e m„ uei G (4+ „)m e m e _ 
i=l r "+i r " + i i r l,2 ^ 

- ELi G4 ";f ' (! + Qe-») 2 - + ae-y ), n, a ~ R (13) 

i=l r i2 ' ' 1,2 J1 

We consider the potential for the case of compactification on an n-dimensional 
torus and assume that the space-time is (4 + n)-dimcnsional, where the n extra 
dimensions are compactified on circles with radius R. Then the (4 + n)-dimensional 
Newton's constant is expressed as^ 

tr(4 +rl ) = ! (, 14 J 

n + 1 i 

where v n = (2nR) n is the volume of n torus and s n = 2ir^~ /T( I1 ±- ) is its surface 
area (T(n) = J™ e" 1 !"" 1 ^). When r ~ i?, for the shape of extra dimensions 
mentioned above, the familiar Yukawa potential is adopted a = 8n/3 and A = R. 
One can notice that the integral diverges as ri 2 — > 0. For n > 2, considering the 
atomic nucleus is not pointlike, we introduce a safe cutoff value r m with the atomic 
nucleus size for the lower limit of the integral. A convenient expression of r m is 
T m = ToA~5 , where A is the mass number of the atomic nucleus and ro is of size 
- 10~ 13 cm. 

If we consider the n-dimensional sphere of radius R as our compactification 
manifold, the Newton constant is 

G(4+ n ) = 2R n G 4 . (15) 

For the case of large extra dimensions, let us compactify the n dimensions y a by 
making the periodic identification y a ~ y a +L. Following exactly the same procedure 
as Ref. |31 one can obtain 

G (4+") = > (. lb J 

^(3+n) 

where V n = L n = (2nR) n and Sd — 27r^/r(-j) is the surface area of the unit 
sphere in D spatial dimensions. From Eqs. (fl4l [T51 IT5|) . since ^p- ~ 2R n ~ g 4lrV ° , 
we only take example for an n-dimensional torus. 

We apply the above wave function in Eq. (JT]) to two-electron atomic and ionic 
systems. In Table 1, we quote the values of the ground state energy for a helium 
atom and its isoelectronic se que nce. For the helium atom we also quote the values 
of Hartree-Fock (HF) theory H^l, the 3-parameter Charatzoulas-Knowleswave (CK) 



function 1221, the variational perturbation results of Pan and the 1078-parameter 
Pekeris wave function For other ions, we compare our calculated values with 
the variational perturbation results of Pan — and Aashamar^^SI. The relative errors 
compared with Pekeris and Aashamar are also given in the last column. 

Table 1. Nonrelativistic ground energies of a helium atom and helium-like ions 
calculated with ip in Eq. {3 in atomic units and their comparisons with other 



I 



I 
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references. For the helium atom, the relative error is compared with the value of 
Pekcris. For the helium-like ions, the relative errors are compared with that of 
Aashamar. 



System 


Wave function 


Ground energy 


Error(%) 


He 


HpHi 
CK m 

Pan EH 
PekerisE21 


-2.90313 
-2.86168 
-2.89007 
-2.89122 
-2.90372 


0.0203 


Li+ 


Aashamar ESI 


-7.27906 
—7 26820 
-7.27991 


0.0117 


Be 2 + 


</> 
Pan 
Aashamar 


-13.65440 
-13.64416 
-13.65557 


0.0086 


B 3+ 


Pan 
Aashamar 


-22.02949 
-22.01973 
-22.03097 


0.0067 


C 4+ 


$ 
Pan 
Aashamar 


-32.40438 
-32.39511 
-32.40625 


0.0058 


N 5 + 


4, 

Pan 
Aashamar 


-44.77929 
-44.77035 
-44.78145 


0.0048 


6+ 


i> 

Pan 
Aashamar 


-59.15416 
-59.14554 
-59.15660 


0.0041 



Furthermore, by making use of the above wave function, we calculate the fre- 
quencies of the spectrum corresponding to the ADD gravity correction. For a he- 
lium atom, we get Av ~ 10 _6 Hz with n = 2, and Av ~ 10~ n Hz with n = 3. 
Although these corrections are much larger than the one calculated from the ex- 
act ISL (At/ ~ 10~ 23 ) for a helium atom, they are still too small to be observed 
in current experiments. So we must find some "large" corrections that cannot be 
simply explained by the inconsistency of the accurate experimental data and Stan- 
dard Model theoretical values. We can choose high-Z helium-like systems (such as 
Ca 18+ and Pb 80+ ) to do the similar calculation. In Table 2, the frequencies of 
the spectrum for two-electron atom and ions corresponding to the ADD gravity 
corrections are presented in detail. However, we should point out that for Pb 80+ , 
since the corresponding Bohr's radius is smaller than that of a helium atom, the 
effect of gravity from the nucleus is much more important. Anyway, if we use the 
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simple perturbation theory, the main correction still comes from the Schrodinger 
term. Because the masses of Ca and Pb are representatively about 40 and 207 times 
of the mass of proton, the corrections of the frequencies for Ca 18+ , Pb 80+ are much 
larger. For Ca 18+ , we have Av ~ 10 _2 Hz with n = 2, and Av ~ 10 _7 Hz with 
n = 3; For Pb 80 +, Av ~ 10°Hz with n = 2 and Av ~ 10~ 5 Hz with n = 3, while the 
corresponding corrections in hydrogen atom are 10 _8 Hz and 10 _13 Hz respectively. 
Also note that the exact ISL, that is, the n = case gives a correction as small as 
Av ~ 10~ 24 Hz for Hydrogen atom and Av ~ 10~ 23 Hz for a helium atom. 

Table 2. The frequencies of the spectrum for two-electron atom and ions cor- 
responding to the ADD gravity corrections 



System 


Az^(Hz) 
(n = 2) 


Ai/(Hz) 
(n = 3) 


He 


-1.1 x 10~ 6 


-1.5 x 10" 11 


Li+ 


-6.7 x 10~ 6 


-9.9 x 10" n 


Be 2 + 


-2.1 x 10" 5 


-3.2 x 10" 10 


B 3 + 


-5.0 x 10~ 5 


-7.9 x 10" 10 


C 4+ 


-9.5 x 10~ 5 


-1.5 x 10~ 9 


N 5 + 


-1.8 x 10~ 4 


-2.9 x 10~ 9 


6+ 


-3.0 x 10" 4 


-5.0 x 10" 9 


F 7 + 


-4.9 x 10~ 4 


-8.5 x 10~ 9 


Ca 18+ 


-1.1 x 10~ 2 


-2.0 x 10" 7 


pb so+ 


-3.1 x 10° 


-7.4 x 10" 5 



We have referred to many references I 24 | 25 | 26 | 27 | 28 | 29 | 30 | 31 | a bout experiment 
values of fine structure interval in helium or helium-like ions, which collected in table 
3. Comparing these experiment values to Table 2, we find that even for Pb 80+ , the 
large shift can not yet been detected because of the precision of the experiment. So, 
the ADD model can not yet be ruled out in this way. We expect some experiments 
will be done in this way in the future. 

To summarize, we have formulated a calculational scheme for detecting the extra 
dimension by using the helium and helium-like systems. To perform the variational 
calculation we have used the variational forms with a flexible parameter k for two- 
electron correlated wave function taking into account the motion of the nucleus. The 
results obtained by us for the ground state is quite accurate. Furthermore, based 
on the results, we illustrate the method numerically by calculating the magnitude 
of the corresponding frequencies. For Pb 80+ , the frequency derived from the ADD 
gravity correction is Av ~ 10°Hz with n — 2 and Av ~ 10 _5 Hz with n = 3. 
The corrections might be used to indirectly detect the deviation of ISL down to 
nanometer scale and to explore the possibility of two or three extra dimensions in 
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ADD's model, while current direct gravity tests cannot break through micron scale 
and go beyond two extra dimensions scenario. 

Table 3. Experimental fine structure intervals for He, Li + , Be 2+ , B 3+ , F 7+ . 



System 


Interval 


Experiment 


Reference 




Vol 


29616943.01(17) kHz 


Ref.Hl 




1>12 


2291161.13(30) kHz 




He 


voi 


29616951.66 T 0.70 kHz 


Ref.[25] 




Vl2 


2291175.59 T 0.51 kHz 






Vol 


29616950.9 T 0.9 kHz 


Rcf.d)] 




voi 


155704.2161(30) MHz 


Ref.HT] 




Vl2 


-62678.3382(27) MHz 




Li+ 


voi 


155704.27(66) MHz 


Ref.Ha 




Vl2 


-62678.41(66) MHz 






voi 


11.5576605(7) cm" 1 


Rcf.[27] 




V12 


-14.892209(1) cm- 1 




Be 2 + 


voi 


11.5586(5) cm" 1 


Rcf.CH 




Vl2 


-14.8950(4) cm- 1 






voi 


16.197573(2) cm" 1 


Rcf.[27] 




Vl2 


-52.661199(4) cm" 1 




B 3+ 


voi 


16.203(18) cm- 1 


Ref.[30] 




Vl2 


-52.660(16) cm- 1 






voi 


-151.2466(1) cm- 1 


Rcf.[27] 




V12 


-957.8487(2) cm" 1 






Vl2 


-957.883(19) cm' 1 


Ref.El] 
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